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Abstract. 

This paper considers work extraction from a quantum system to a work storage 
system (or weight) following reference [1]. An alternative approach is here developed 
that relies on the comparison of subspace dimensions without a need to introduce 
thermo-majorisation used previously. Optimal single shot work for processes where a 
weight transfers from (a) a single energy level to another single energy level is then 
re-derived. In addition we discuss the final state of the system after work extraction 
and show that the system typically ends in its thermal state, while there are cases 
where the system is only close to it. The work of formation in the single level transfer 
setting [1] is also re-derived. The approach presented now allows the extension of 
the single shot work concept to work extraction (b) involving multiple final levels 
of the weight. A key conclusion here is that the single shot work for case (a) is 
appropriate only when a resonance of a particular energy is required. When wishing 
to identify “work extraction” with finding the weight in a specific available energy 
or any higher energy a broadening of the single shot work concept is required. As a 
final contribution we consider transformations of the system that (c) result in general 
weight state transfers. Introducing a transfer-quantity allows us to formulate minimum 
requirements for transformations to be at all possible in a thermodynamic framework. 
We show that choosing the free energy difference of the weight as the transfer-quantity 
one recovers various single shot results including single level transitions (a), multiple 
final level transitions (b), and recent results on restricted sets of multi-level to multi¬ 
level weight transfers. 
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1. Introduction 

The neat characterisation of general classical non-eqnilibrium processes in terms of 
fluctnation relations la El SI m has rapidly advanced the general understanding of 
thermodynamic processes and properties at the mesoscopic scale. Work, in particular, is 
a thermodynamic quantity of interest and the stochastic fluctuations of work done on a 
system are captured in the Jarzynski relation [3]. The fluctuation relation approach 
has been extended to quantum systems where probabilistic energy transfers of the 
system that undergoes unitary evolution are associated with the fluctuating work done 
on the system [HiEiEiiniiininiiis]. Again the work for this quantum scenario can be 
captured in a quantum Jarzynski relation. On the other hand, thermodynamic processes 
for a quantum system can be studied in a setting where the system interacts with a 
heat bath and a work storage system (or weight) undergoing global unitary dynamics 
na m m Eg EHi EHi [lEi tta E9] .' Operation on the system then results in a change of the 
work storage system’s state and it is that change that is here associated with “work”. 
These approaches are referred to as “thermodynamic resource theory” and “single shot 
thermodynamics”. Recent papers, e.g. PESIES], derive upper bounds on the amount 
of work that can be drawn from a quantum system that starts in a non-equilibrium 
state in a “single shot”. 


a) 


w 


(c) 




Figure 1. Sketch of the transition of the weight (a) from a single energy eigenstate 
(of energy 0) to another single energy eigenstate (of energy w), (b) from a single energy 
eigenstate of energy 0 to a set of energy eigenstates of energy [iD, re -I- J], and (c) from 
a range of energy eigenstates to another range of energy eigenstates. 


The single shot work done by the system is here associated [T] with the transition 
of the weight from a single energy eigenstate (of energy 0) to another single energy 
eigenstate (of energy w). This situation is sketched in Fig. The proof of the 
bounds provided in [1] relies on established mathematical concepts from quantum 
information theory and a new majorisation concept called “thermo-majorisation”. 
Researchers, in particular those outside of quantum information theory, may find the 
proof mathematically heavy and are unable to follow the detailed logic. This technical 
difficulty overshadows the interpretation of the results and may hinder their further 
development and formulation of experimental tests of single shot results. 

Here we aim to re-derive single shot work extraction limits while keeping the 
technical side as simple as possible. The hope is that stripping the discussion from 
some of the jargon will allow to focus on the physical meaning of the results, clarify 
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the situation they describe, and develop the argument further to adapt to different 
physically relevant scenarios. 

The paper is organised as follows. The known work extraction result of [1] for 
scenario (a) in Fig. is re-stated in Section]^ and then re-derived in Section]^ The hnal 
state of the system after this work extraction procedure is explored in Section Also 
re-derived, in Section is the work of formation which was hrst identihed as different 
from the extractable work in [T]. The extension of single shot work extraction to general 
transfer processes is discussed thereafter. Section is concerned with the extractable 
work allowing (b) the transition of the weight from a single energy eigenstate of energy 
0 to a set of energy eigenstates of energy [w,w -f 5]. Section identihes a transfer- 
quantity that characterises (c) the transition of the weight from an arbitrary range of 
energy eigenstates to another range of energy eigenstates. Finally, the hndings and open 
questions are discussed in Section [8j 

2. Known resnlts on maximal single shot work extraction 

A key result in single-shot thermodynamics is the identihcation of a maximal work uma 

wT- = Fr{ps)-F{Ts). ( 1 ) 

that can be extracted with success probability 1 — e from a system starting in a state 
ps under so-called thermal operations flj. Here F{ts) := —^ InZ^ is the standard free 

energy associated with the thermal state ts := ^ system Hamiltonian Hs at 

inverse temperature (3, with Zs = tr[e“^^®] the partition function. (In pQ the system is a 
qubit with the excited state energy tuned such as to have exactly the optimal work value 
that can be gained, i.e. Hg = 0 |0)5(0| -f |l)s(l|.) The other quantity, 
is a generalised free energy applicable for the non-equilibrium state ps [I] which will 
be detailed below. This maximal work is valid for initial states ps that are diagonal in 
the energy basis and it is at least a lower bound on the maximal work for non-diagonal 
states mmm- 

This result is derived within the thermodynamic resource theory setting [ISlElttl 
[la nn [ig ng [m [19] involving three components: the system of interest, S, a bath H, 
and a weight (or work storage system) W. In the simplest case, the Hamiltonian at the 
start and the end of the process is assumed to be the same and the sum of the three local 
terms, H = Hs + Hb + H^r. The weight is assumed to have no degeneracies, i.e. all its 
energy eigenstates have different energies Ew- The system’s degeneracy is not restricted 
and we denote the multiplicity of each energy Eg by Ms{Es) and label each of them by 
gs{Es) = 1,..., Ms{Es). The degenerate bath levels at energy Eb are labeled by /b^Eb) 
and their multiplicity is assumed to be exponentially growing with inverse temperature 
Jb^Eb) = 1,Mb{Eb) = 1,...,Mb( 0) exp(/3F^s). This exponential growth of the 
degree of degeneracy (or the density of states) generically results for all systems made 
of many similar subsystems that feature short range interactions izu. We note that it 
only holds for a possibly large but hnite energy regime. 

Thermal operations have been dehned [I] as those transformations of the system 
that can be generated by a global unitary, V, that acts on system, bath and work 
storage system initially in a product state ps Z) tb Z) \E^^)\y{E^w\, with the bath in a 
thermal state, tb = ^~^~i and the work storage system in one of its energy eigenstates. 
Conceptually, global unitaries V describe the operation of a “work extraction 
machine” that aims to extract as much energy as possible to the work storage system. 
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Perfect energy conservation is imposed by requiring that the unitary may only induce 
transitions within an energy shell of total energy E = + Es + Eb- This is equivalent 

to requiring that V must commute with the sum of the three local Hamiltonians. 

The generalised free energy of the non-equilibrium state ps is defined as 

Fr(Ps) := In ^ (2) 

^s^gs 

Here h{Es,gs,^) is a binary function that determines whether a particular energy 
eigenstate, \Es,gs), is included in the summation or not. The value of h depends 
on ps and on the failure rate e that is being accepted for the work extraction. The exact 
dependence will be discussed further in Section The optimal work stated in Eq. ([^ 
now corresponds to the following task: For the weight initially entirely in its ground 
state (with energy E^^ = 0) the full system is transferred to a final quantum state such 
that the probability for the weight to be found in an energy eigenstate with energy E^ 
is 1 — e where 1 > e > 0, see Fig. This process is associated with the “lifting” of the 
weight and the energy difference experienced by the weight is identified with “extracted 
work”, w := E^ — E^\ For the special case of perfect work extraction, i.e. e = 0, the 
summation in Eq. (j^ includes all energy eigenstates \Es,gs) that are populated in the 
initial state, i.e. h{Es, gs,^) = 1 when {Es, gs\ps\Es, gs) > 0, and 0 otherwise. The 
maximal work is then 

=-^ lntr[rsflp3], (3) 

where Hp^ is the projector on the support of ps- 

3. Work extraction bounds from limits on probability transfer 

The derivation of the maximal extractable work presented in [T] rests on an analysis 
of state dimensions which is combined with the newly introduced concept of thermo- 
majorisation. This is a variation on majorisation, an important tool in the study of 
doubly stochastic matrices [22l |23l El] and quantum channels in particular [25l 126] . 
Here we present a derivation of the same result without a need of invoking thermo- 
majorisation explicitely. We hope this re-derivation is more straightforward to follow 
for researchers wishing to familiarise themselves with single-shot thermodynamics. 

The rational of the presented approach is to compare dimensions of the involved 
projective subspaces to conclude what transformations on the system are possible in 
the setting given and what maximal work extraction they enable. The re-analysis opens 
avenues of generalising the previous result, valid for transitions of the weight from single 
energy level to single energy level, to processes where the weight transfers between 
multiple energy levels, which we will explore in Sections Eland The work required to 
form a non-equilibrium state is also described, in Section]^ 

The desired thermal transformation is of the form 

g ■= PSB ® |0)w(0| (TSB ® \w)w{w\ =■ v'{y), (4) 

where the system and bath start in a fixed state psB and end in some state asB, while 
the work storage system lifts from a single energy level, |0), to another single energy 
level, |tc), by an energy tc > 0, see Fig. Such a transformation, however, turns out 
to be impossible for a large class of initial states, psB, namely those which have full 
rank. Thus one allows transformations under which the l.h.s. of (|^ is transformed into 
the r.h.s. only up to a success probability 1 — e, as already introduced in Sect. 
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Since the global unitary, V, is fully energy conserving we may treat the global 
dynamics for each global energy shell E separately and later combine their contributions. 
Note that those energy shells are projective subspaces and so are the individual energy 
levels that make up the shell. We denote the local projectors dehned by their associated 
local energies Es,Eb and Ew by IIs{Es),IIb{Eb) and IlwiEw)- 

To hnd the maximal work we hrst consider a probability 1 which falls 

initially into a projective subspace 11;^^ of dimension d®-, i.e. P,® = tr[r 7 n;^;] for the 
initial state rj. The key task here is to decide if this probability can be transferred 
entirely into a projective subspace under unitary transformations V. This is of 
course possible if 

tr|nf„ 1/ nf., r, nf„i v'] i tr|nf„, r, ngj = pr (5) 

holds. Rewriting this condition using the non-trivial eigenstates |n) of R r/11^; Rf, 
i.e. RlfR r^IfR V^\n) = PnW) with > 0 for n = 1,..., dR, implies 

dL ^ ^ dg, 

(n|nfjn) = = Pifi (6) 

n=l n=l 

and therefore one must have {nllflj^ln) = 1 for n = 1,..., d®,. This is only possible if the 
dimension of the space into which the probability is mapped, df^^ := rk[rifj^], is at least 
the same as the initial dimension, i.e. the condition becomes 

dfn > dfni- (7) 

If not the whole probability PR but a slightly reduced probability (1 — e)PR is to be 
transferred, then one may replace flR in Eq. (|^ with any projector nj^i(e) that fulhls 

( 8 ) 

Dehning the smallest initial subspace dimension as dfffe^rj) := miupB / dim(nP(e)) < 
dR, the hnal space dimension condition, when allowing a small error probability e in 
the transfer, relaxes to 

din>dfr,i{Fri)- (9) 

This constraint is now the key condition to establish what the maximum extractable 
work is for the thermal operations specihed in Eq. Q. 

Since R cannot mix energy shells, it is necessary to check the condition for each 
energy shell E. We interpret the projector used above as the operator that projects 
on the non-trivially populated subspace of the hnal global state, ass ® with 

energy E. By construction one has the relation 

n|n = ^ ^s{Es) ® ^b{E — Es — w) ® tlw{w). (10) 

Es 

We can now identify the dimension of n|(,, 

dl = Ms{Es) Mb{E-Es-w) = J2 MsiEs) Mb{E) (11) 

Es Es 

= MB(E)e-^'‘Zs, 


( 12 ) 
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where Zs = 'YIe ^s{Es) is the thermal equilibrium partition function at inverse 

temperature [3 for the system. The second step is to consider the initial subspace 
projectors in the energy shell E, 

flE = MEs) ® Ub{E - Es) ® Ilw(O). (13) 

Es 

In the following we aim to specify the dimension dfffe^rj) of a sub-projector of flj^j in 
whose associated subspace lives the fraction 1 — e of the initial state {rf) population (PifJ 
in energy shell E. For e = 0 and psB a full rank state the dimension of itself is just 

4i(0) = Mb{E - Es) = Mb{E) Zs. (14) 

Es 

However, to determine the initial dimension d^fe, rj) for a finite failure rate e > 0 we will 
need to specify the global initial state rj a little more. While not the only choice, here 
we consider the class of initial global product states previously discussed [BIISIIISIE], 

ri = ps ® tb ^ \0)w{0\ (15) 

where ps is a system state that is diagonal in the basis of the system Hamiltonian and 
the bath is in a thermal state tb- We will label the eigenvalues of ps when diagonalised 
in the system energy eigenbasis by X{Es, gs) where Es is the corresponding energy and 
Qs the degeneracy index. By construction all non-zero eigenvalues of p in the energy 
shell E, i.e. the non-zero eigenvalues of 11 ®, p H®,, are then 

I3(E-Es) 

f'^{Es, gs, fs) = YEs, gs) -^-• (16) 

Here fB = 1, Mb is the degeneracy index of the bath at energy Eb = E — Es] but 
it does not affect the magnitude of r/’s eigenvalues. The eigenvalues thus have a high 
multiplicity, given by the bath multiplicity Mb- The eigenvalues can be re-labeled with 
a = 1 ,..., d®;( 0 ) as r^{Es, gs, fs) = arranged in decreasing order, 

rf > rlf > r®..., see Fig. The spectra of different energy shells E, as sketched in 
Fig. differ only by a factor of e“^® for the individual values while their multiplicity 
Mb differs by a factor of e^®. The population probability of the global state p in energy 

shell E is given by = ^ifi- 

Now we are ready to identify the dimension dffe,p) of the subspace within the 
energy shell E where most of the initial probability, P®, (1 — e), lives. In the sum above 
one simply has to add up the a to dE{e,p) < d®i( 0 ) which must be chosen such that 
the probability reduces to the desired level, 

Y1 rf > p,5 (1 - t). (17) 

CK=1 

Comparing with Eq. 0 the final dimension is given by 

v) = ^ Mb{E- Es) h{Es, gs, e) (18) 

i^s,gs 

where ^(P^, ( 75 , e) determines which terms are included in the sum, as visualised and 
explained in Fig.l^ An h-value of 1 (0) indicates that a whole block of global eigenvalues 
all corresponding to system state \Es,gs) is (not) included in the summation, while a 
fractional h-value indicates that a fraction of the eigenvalues in a block is included, see 

Fig. II 
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Figure 2. Eigenvalues of global initial state rj in energy shell E given in Eq. (16) 
arranged in decreasing order by their index a. For any particular system state \Es,gs) 
there is a whole block of global state eigenvalues of the same magnitude that arises 
due to the bath’s degeneracy Mb{E — Es) at energy Eb = E — Es- The full initial 
probability, PF, is indicated by the red box. The slightly reduced probability that is 
to be transferred, PF (l — is indicated by the green box. To make up the slightly 
reduced probability blocks of probabilities are either fully included [h = 1), fractionally 
included (e.g. h = 3/5), or not included {h = 0) in the summations Eq. (18) and 
Eq. (pi). 


The value of df^fe^rj) must be chosen as the smallest integer that fulhlls (17). It 
can happen that e has a value that would require to split up an individual eigenvalue 


in order to fulhll (17) as an equality, i.e. + aC+i = PF (1 — e) where n is 


an integer and 0 < a < 1 a real fraction. In this case the eigenvalue must be fully 
included, i.e. h must be chosen as the larger proper fraction available. For example, 
having to split the 3rd of 5 eigenvalues in a block {Es,gs) the associated h{Es, gs,^) 
must be chosen 3/5 rather than 2/5, in order to guarantee that the failure probability 
is strictly less than e. 

For a map that is identical to (|^ up to a failure probability e we have obtained the 
initial and hnal subspace dimensions for each energy shell E. They allow us to check 
the maximal work extraction condition Eq. (1^ for each E. One has 


dl = Mb{E) e-f^^ Zs >Y. ^b{E) h{Es, gs, e) = dgfe, g) (19) 

Es>9S 

w < - ^ In ^ h{Es,gs,e) - E{t) = 




The last line shows that because of the trivial dependence of d^fe,g) on E, the latter 
drops out and the expression for the maximal extractable work stated in Eq. ([^ is 
recovered [1]. Note, that the dehnition of h has been extended from a binary function, 
used in Eq. ([^, to include rational fractions. The optimal work described by Eq. ( Ip ) 
is thus marginally larger than the work described by Eq. To recap, this work is 
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maximal under variation of the global unitaries V that realise the map (|^. The bound 
is tight when e is chosen such that no splitting of eigenvalues would be necessary. 

Note, that using the eigenvalues of the thermal state, ts, denoted by t(Es) := 
the maximal extractable work ([^ can also be written as 

= -4 In E t{Es)HEs,gs,0- (20) 

-EsiSS 

Clearly, if h was 1 for all its arguments, i.e. all the eigenvalues in Fig. have to be 
included, then the maximal extractable work is 0. For e = 0 no probability can be 
lost and h = 1 for all energy levels with non-zero population. Work can thus only be 
extracted if the initial state has rank less than the system Hilbertspace dimension or 
e > 0 or both. 


4. The final system state 


The work extraction protocol, Eq. Q, has mapped r] —)• rf'iV). It is natural to ask in 
what reduced state as = tr^vvlh^^)] fhe system is left as a result. We will here answer 
this question. 

A crucial point to note is that not a single, unique unitary V corresponds to 
maximum work extraction but that there are inhnitely many. Any unitary that maps 
the projective energy subspaces n®j(e), with dimensions df^fe^rj) calculated in Section 
[i, onto the projective subspaces ffl'j^ will give rise to maximum work extraction. The 
large set of unitaries with this property forms a subset of all allowed unitaries, i.e. those 
that commute with H. As a hrst step we look at the average hnal state, obtained 
through the application of each of these unitaries and integration over the respective 
Haar measure, which we denote by ((■)). Details on integration over Haar unitaries can 
be found e.g. in [27]. Note that even though the unitaries V commute with if, 
itself need not be diagonal w.r.t. to the product basis of local energy eigenstates that 
dehnes the energy shells. While H and rj^V) share a simultaneous eigenbasis there is no 
guarantee, due to H having degeneracies, that this is the product basis of local energy 
eigenstates and thus rj'iV) can have off-diagonal elements in the degenerate blocks of if. 
For the diagonal part of r]'{y) w.r.t. the above basis consisting of products of energy 
eigenstates of individual subsystems, '? 7 diag(D), one hnds 

((hdiag)) « ® 0 \w)w{w\, (21) 


Appendix Ajfor details of the derivation. This is a factorized state where the system 


see 

and bath are each in thermal states at inverse temperature (3. However, this result in 
itself is not sufficient to conclude that typically an almost thermal state of the system 
results from an individual V. One needs to calculate also the relative variances of 
Es, gs)- The latter are the portions of the diagonal elements at Es,gs of the 
hnal reduced state of the system as that correspond to the parts of the total hnal 
state living in energy shell E, i.e., (Es, gsWslEs, gs) = Y^e P&n{E, Es, gs)- Hence, if 
these relative huctuations are small one will typically get diagonal elements of as that 


resemble their averages, which in turn coincide, according to Eq. (21), with those of a 
thermal state. Concretely we obtain for those relative huctuations, see [Appendix A 

AFfin(E, Es, gs) ^ _1_ 

{{PM(E,Es,g3))} ~ s/Mb(E - Es-w)' 


( 22 ) 
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which is indeed very small due to the large multiplicity of the bath Mb{E — Es — w) = 
Mb{E) for energies E » Es + w. Off-diagonal elements are considered 

separately, as detailed in Appendix A but it turns out that the variance to their 

Thus, for a unitary drawn at random 


close to the average 


vanishing mean also scales as -. 

yjMB(E-Es-w) 

all elements of a typical final system state as are , - 

Mb(E-Es-w) 

thermal state if the bath is large. The entire above reasoning is along the lines of what 
has become known as “quantum typicality” |2Hl |29l [30l |3T] . 

Having realised that optimal work extraction can result in a multitude of final 
states, with typical state being close to the thermal state, one may wonder why it is 
not possible to extract more work in those rare instances where the final state happens 
to be fairly different from the thermal state. Here it is important to note that the final 
state need not to be factorized between system and bath (c.f. Eq. (21) only describes 
the diagonal part). Indeed, it is known [32l |33l |3ll |35l [36] that only a negligible set of 
the final states will be factorized, especially when the degeneracies of the bath B are 
large. However, to repeat the same work extraction process requires an initially fully 
factorized state, thus preventing repetition and further work extraction along the same 
route for the majority cases. 

Finally, a concern may be that one could introduce a “fresh bath” that is 
uncorrelated with the final system in order to extract even more work from a non-thermal 
hnal system state. Here it is important to realise that the work extraction process with 
the first bath left the system not necessarily in a thermal state, but certainly in a state 
with full rank, as- This is due to the fact that maximum work extraction requires a non¬ 
zero population of every dimension spanned by n|^, see Eq. (10). As can be seen from 
Eq. (|^ a full rank system state means that there is no potential for further “perfect” 
work extraction, i.e. moving weight population from |0)^y —)■ \w)w with certainty 
(e = 0), even if new baths are brought in. If the first work extraction process is imperfect 
(e 7 ^ 0), i.e. the weight ends with population in |0)iy and \w)wi in general weight and 
system will be correlated, thus also violating the initial factorisation condition, see (|^, 
for a second work extraction process. Nevertheless, for completeness, suppose that the 
result factorises and the weight is chosen such that it has additional energy levels [^ 2 ) 
and |ta -|- 1 ^ 2 ), where W 2 is adjusted to match W 2 = wf'^^^as) of Eq. (20) applied to as- 
One can now indeed perform a second imperfect work extraction which may very well 
increase the mean energy in the weight system. However, this will lead to a weight level 
population of at least three of the four eigenlevels, \Q)w)\w)wA'^‘i)wA'^Ew 2 )w- This 
scenario is outside of what is allowed in case (a), see Fig. [^, but will be discussed as 
case (b) in Section If, on the other hand the first work extraction process is a perfect 
one, weight and system necessarily factorise and a second, imperfect, work extraction 
can be added. This will result in overall imperfect transfer of |0) to as a 

result of the total process instead of raising perfectly to in the first step 

alone. So one could have obtained the same result if one had allowed for imperfect work 
extraction in the first place since an increase of e always allows for an increase of work, 
y^max(^^) ^ This sliows that there is no second law conflict with a final state 

being non-thermal. 


5. Minimal work cost of formation 

A key finding of the single shot thermodynamics approach presented in [T] is the 
realisation that, in analogy to entanglement of formation and distillation, the work 
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required to form a non-equilibrium state may differ from the work that can be extracted 
from that state. Here we will now re-derive the minimal cost of formation, w™™, of 
engineering with certainty (e = 0 ) a diagonal state as from a thermal state ts = ^ , 
under thermal operations. 

Considering again global unitaries V that commute with the Hamiltonian, the 
desired operation is 


r] : = 


Ts <S)Tb® \w)w{w 



(Tsb ® |0)vi/(0| := 


(23) 


with any final state permitted such that trB)^^^] = as- Using the same reasoning that 
lead to Eq. (16) the eigenvalues r^{Es, gs, fs, Ew) of the global initial state in a 
particular energy shell E are given by 


^-g(E-w) 

r^{Es,gs,fB,Ew) = „ „— Sew,w, (24) 

where gs and fs are the degeneracy indices labelling the various energy eigenstates for 
Es and Eb = E — Es — Ew, respectively. By construction the global eigenvalnes are 
only non-zero when the weight has energy Ew = w and the non-zero eigenvalues in shell 
E are all equal. 

The eigenvalues of the desired final system state <75 which is assumed diagonal 
in the system’s energy eigenbasis {|E 5 , 5 f 5 )} of Hs are denoted by s{Es,gs) = 
t>T^[<TsB \Es, gs){Es, gsW- The contributions stemming from the energy shell E to the 
diagonal elements of the final system state as rnay thus be written as 


s^{Es,gs) = - EsjB\{Es,gs\v' \Es,gs)\E - EsJb)\0), (25) 

fB{E-Es) 

where the snm is performed over all bath degeneracy indices at energy Eb = E — Es- We 
note that since the initial state has no off-diagonal elements between different system 
energies, say, Es and Eg, and the nnitary V commutes with H no such off-diagonal 
elements can be generated under V. There may be off-diagonal elements within a single 
energy snbspace Es, i.e. {Es, gsWs \Es, g's) 7 ^ 0- These off-diagonal elements can always 
be “avoided” though by choosing the basis {IE 5 , gs)} for the system subspace of energy 
Es appropriately. Thus, w.l.o.g. the state as can be assumed diagonal with diagonal 
elements as given in (25). 

Since V is unitary, not only the eigenvalnes of the initial state rj (in energy shell E) 
bnt also those of the final state g' (in energy shell E) are given by Eq. (24), i.e. they are 

either 0 or ^ Zg ' summands on the r.h.s. of (25) are just expectation values of 
the state g' and since any expectation value of a state is upper bounded by the largest 
eigenvalue of the respective state we find from (25) 

Q-h{E-w) 


!>‘^(E3,9s)< Y. 

fsiE—Es) 


ZsZ 


B 


^-h{E-w) 

Zs Zb 


Mb{E - Es) = 




(26) 


-‘B 


The probability of the final state as to be found in state \Es, gs), denoted by s{Es, gs), 
is then obtained by summing over contributions from all the energy shells. 


s{Es,gs) = '^s^{Es,gs) < 

E 


Q-hiEs-w) 

Ts 


(27) 


where MB{E)e = Zb was used. Importantly this ineqnality must be fulfilled for 
all pairs Eg, gs- Using the eigenvalnes of the thermal state of the system and rearranging 
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one obtains 


wz? ^ ^ 1 


(28) 


as the condition on the work cost of formation. One can see now that the energy 
eigenstate \E*g,gg) with the largest ratio ^ F is the one 


that constraints the whole transformation. The tight lower bonnd on the work cost of 
formation is then 


w > — In /i =: 
P 


(29) 


fi can also be expressed as /i = min{A : s{Es,gs) < l^t{Es)'^Es, gs} or 

p = min{A : < A r^} (30) 

for as the desired diagonal state of the system. The work cost Eq. (29) for the ideal 
map Eq. (23) can easily be extended to allow hnal states in an e vicinity of the desired 
state as- /i is then replaced by = info-| min{A : a^ < Xts} where a^ is close to as, 
||(t| — (T 5 II < e, where || • || is the trace norm m- 


6. Generalisation of work extraction to multiple levels in the final state of 
the work storage system 

We have seen in the previons section that the single shot concept gives sensible results 
for work extraction when single levels of the weight are being considered (case (a) in Fig. 
1). While for truly microsopic systems this may be a suitable analysis, case (a) is not the 
right picture for mesoscopic or macroscopic work extraction process. It is not feasible to 
lift a work storage system from a true quantum energy eigenstate to another true energy 
eigenstates. The density of states is so large that it will be practically impossible to 
pick a single eigenstate of a macroscopic system This indicates that an extension 
of the single shot thermodynamic framework to transitions between multiple levels, see 
Fig.[T]3 and [^, is desirable. 

The approach presented above now allows a straightforward generalization to 
thermal processes where the work storage system is transferred from the ground state 
to not just a single energy level, but an interval S of energy levels, see Fig. [pz. To do 
so we will now turn to a harmonic oscillator work storage system with energy spacing 
AE, similar to the one considered in mi. and differing from the qubit work storage 
system used in [1]. For simplicity we will assume that also the system and the bath 
have oscillator-like equidistant energy eigenlevels of integer multiples of AE. Again, 
the weight shall have no degeneracies whatsoever while the system and bath again have 
degeneracies Ms{Es) and Mb{Eb) at their energies Es and Eb, respectively. 

Work extraction of an amount of work is now identihed with the transition (with 
probability 1 — e) of the work storage system from its ground state to any state living in 
the energy subspace [we, M)e + 5], see Fig. 03- The ideal transfer (e = 0) is now associated 
with the map 

V'■= Ps®'Tb®\^)w{^\ -P (^SBW (31) 

where the hnal state of the weight aw := tr sbWsbw] is now mixed with support in the 
interval specihed. Here the general case of entangled states between weight and the rest 
is allowed. If the role of these correlations is of interest, then restrictions to product 
states could be considered. 
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To account for the different (larger) space of possible final states one simply has to 
restart the calculation of nf„ and as detailed in Sec. and enlarge the final Hilbert 
space in order to comprise all final global states with the weight in the specihed local 
energy subspace. This “new” final Hilbert space can be described as a sum over many 


“old” ones. Restarting from (11) one can introduce a pertinent summation over weight 
energy levels so that the new reads: 


w+5 


^-0Ew 


(32) 


dt = Ms(E) Zs Y. e 

Ew=w 

The sum on the r.h.s. now depends on the level spacing AE which did not play any 
role so far, 


W-\-S 

Ew=w 


(w+S) / AE 


^—fSEw _ 


E 

n=w/AE 


—/3AEn _ 


_ g-/3((5+AE)^ 
1 - ^ 


(33) 


where n is a natural number that labels the eigenvalues in the allowed hnal interval 
Substituting this for e~^^ in ( pT| then leads to 

1 1 _ g-/3(5+Ai<;) 


^max _ ^max 




1 - 


(34) 


Increasing the width 5 of the allowed final energy regime for the weight beyond a single 
energy level {6 > 0) leads to an additional term with respect to the single final energy 
level situation characterised by Eq. ([^. For large energy intervals (d) and small energy 
spacings {AE) the last expression quickly converges to 

(35) 


^max _ ^max 




Expression (35) holds under the conditions that compared to the “thermal energy”, | 
the final energy range is large, d >> ^ and the energy spacing of the work storage 
system, AE, is small AE « T Physically this means that the final energy range d 
has to comprise a large number of energy eigenstates of the work storage system. 

We note that the extractable work for case (b) in Fig. is increased in 

comparison to the extractable work of case (a) in the limit considered {AE « 
1//3). The surplus on the r.h.s. does not depend on d as long as the above conditions are 
fulfilled. However, the work does vary with the level spacing and can go up to infinity 
for AE —0. 

To judge whether these conditions are fulhlled for meso/macro-scopic work storage 
systems, it is instructive to consider realistic numbers. At room temperature the thermal 
energy is ^ ~ 2.5 ■ 10“^ eV; about the energy of a single optical phonon in a solid. 
Thus allowing a final energy range of d for the macroscopic weight system that is large 
compared to the energy of a single phonon, condition d >> ^ is fulhlled. For the work 
storage system one can imagine a pendulum with eigenfrequency 1 Hz which results in 
an energetic level spacing of AE k, 4-10“^® eV<< ^ thus fulhlling the second condition. 
Level spacings of many meso/macro-scopic work storage systems, such as batteries, are 
even smaller by several orders of magnitude. These example numbers indicate that both 
conditions will generally hold for meso/macro-scopic systems. 

Maybe surprisingly the additional term, — ln(/3AF)//3, is an up-shift as sketched 
in Fig. 03, implying a work above that of Eq. ([^. This seems to indicate a conhict of 







From single-shot towards general work extraction in a guantum thermodynamic framework 13 


Eq. (35) with the second law, but this is only a problem in a naive sense. To properly 
discuss second law violations one would have to build a cyclic machine that enables work 
extraction. It is possible to close the desired work extraction transformation, Eq. (31), 
to a cycle, however, it is not possible to run the cycle again. This is because the weight 
started in an energy eigenstate, |0)vi/, while ending in a reduced state aw which is not 
an energy eigenstate. Therefore there is no true second law violation as the process 
cannot be repeated using the same weight. Exploring if there is a link between the 
apparent violation of the second law in Eq. (35) and the existence of weight-system or 
system-bath correlations [33 ESI ESI in the final state in Eq. ( |3l| ) is an interesting future 
avenue. 

Finally, the implication of Eq. (35) is that one can extract more work when the 
final weight state is allowed to live in a range of energy levels (b), rather than when 
restricting to a single level (a). While this is mathematically sound one may find this 
physically counterintuitive as single shot work is colloquially often associated with a 
worst case scenario. The sketch in Fig. intuitively suggests that the worst case result 
for jumping from the ground state to a range of energy states is a single level transition 
to the lowest level, implying the lowest energy in (b) is the same as the single energy 
in (a), w = w. However, this is not correct - the salient point is that populations of 
energetic levels higher than w would not count as “success” in the situation depicted 
in Fig. [^. To achieve a high probability of success, 1 — e, in (a) the population needs 
to be concentrated in just that one level and this is what leads to a lower amount of 
extractable work in comparison to case (b), i.e. w > w. This new result may appear 
to contrast with a recent paper 001 that shows that allowing population in higher 
levels should not change the “predicted” work. However, the work definition used there 
follows traditional statistical physics concepts, without including a weight explicitly and 
involves sequences of shifting energy levels in the system Hamiltonian and coupling the 
system to a bath na M\. The work concept thus differs signihcantly from the work 
concept used here. 

Now the question is how significant it is if the weight ends up exactly with energy, 
say w, or with an energy in a range, say [w, w S]7 The second case is of relevance for 
many mechanical processes in physics and engineering where gaining work is the key aim, 
and obtaining a little more energy would be judged positively. There are other situations 
where resonance with a particular energy value, i.e. case (a), is key for performance, 
for instance, this may be the case in biological energy conversion processes, such as 
photosynthesis [32] • We propose that the work “extracted” in such resonance processes 
may be described as “resonance work” or “matched work”. To decide which approach 
to use to calculate work, e.g. the single shot approach with single level transitions (a) 
or multiple level transitions (b), or other approaches, one has to first identify what is 
really needed - energy of a certain amount or energy above a certain threshold. 


7. Generalisation of work extraction to multiple levels in the initial and 
final state of the work storage system 

We found that the work definition based on a “single energy level to single energy level 
transition” of a weight, case (a) in Fig. is too narrow and at odds with our notion of 
classical work storage systems. A physical interpretation was however possible in the 
context of matching a particular energy, i.e. in a resonance situation. For the physically 
motivated extension to multiple final levels, case (b) in Fig. that corresponds to work 
extraction in the sense of “at least this much or more energy” the predictions came out 
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Figure 3. Weight energy probability distributions, P{Ew), over weight energy, Ew- 
Sketched are (a) (green) a single level to a single level transfer, (b) (orange) a single 
level to a range of levels transfer, (c) (blue) a range of levels to a range of levels transfer, 
cf. Fig. [2 The particular choice of the distributions is not generic - the orange box 
shows the highest entropy state for that particular range, and the initial and final blue 
distributions are chosen the same shape, implying that these distributions have the 
same entropy. 


unphysical. So one is left with the unsatisfactory situation that there is no consistent 
single shot work extraction concept for general initial and final states of the weight. 
Such extension issues have been recognised by others and motivated a number of recent 
publications that provide bounds on the single shot work from majorisation |1D] and 
fluctuation theorems [13], address the shifting of general weight probability distributions 
ini, and employ axiomatic approaches to identify sets of monotone functions that can 
act as work quantihers for specihc situations [Ml Il5] . 

To try and rectify this situation we here proceed to allow multiple energy levels 
in both, the initial and final state, as indicated in Fig. and Fig. &■ Here one 
now faces the problem of defining work on the basis of an initial and a final energy 
distribution of the weight. When the weight transfers from multiple levels to other 
multiple levels without a means to distinguish any of these transitions, there is no 
single energy difference AE^y that can be straightforwardly linked to work. Even if 
one considers the multitude of energy differences associated with the different single 
level transitions, it is challenging to formulate a “worst case scenario” for the reasons 
discussed in the previous section. We will here not attack the complicated problem of 
defining single shot work for transfers between multiple levels of the weight, see Fig. 
and Fig. [^, and finding the optimum work. Instead, our aim here will be to find a 
transfer-guantity, which we denote (w), that indicates whether an initial weight energy 
distribution can or can not be transformed into a hnal distribution. To analyse the 
general multiple level scenario we consider again a global unitary {V) commuting with 
the sum of local Hamiltonians as before and no restrictions on their spectra. V maps 
an initially factorised state to some final possibly correlated state, 

r] ■■= Ps®tb® (Tw , (36) 

where the initial bath state is assumed thermal and the weight’s final reduced state is 
cr^ := AseVi']- Various transfers between different energy distribution of the weight are 
sketched in Fig. 

The average energy change of the weight on its own cannot be signihcant to limit 
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thermodynamic transformations. For example, a weight starting in the ground state 
can have its average energy raised just by coupling it to a heat bath. Intuitively, the 
transfer-quantity should reward energy increase while punishing spreading the energy 
across different energy levels of the weight. It is known that the free energy has exactly 
such a balancing property and this motivates us to dehne the transfer-quantity for 
bringing a weight from state aw to state as the free energy difference of the weight, 

(w) := AFw = AUw - TASw- (37) 


Here AUw = t,T[Hwicr[Y — aw)] and ASw/ks = —ir[a'^r Incr^y] -|- tT[aw In ctw] are the 
average energy change and entropy change of the weight and T is the temperature of the 
initial bath state. Using sub-additivity of the von Neumann entropy it is straightforward 
to show, see Appendix B, that the transition Eg. (36) can only be possible when the 


transfer-guantity obeys at least the following inequality 

{!»> < F(ps) - F(p's). (38) 

Here F{ps) — F{p'g) = —AFs = —AUs -\-TASs is the free energy difference of the two, 
in general, non-equilibrium system states ps and p'^ =:= trswiv'], with AUs and ASs 
dehned analogously to above, using instead the system Hamiltonian Hs- Note, that 
inequality Eq. (38) is necessary but not sufficient, i.e. there are examples where this 
inequality is fulhlled and the transition may still not be possible. Equality in Eq. (38) 
can only occur when the hnal global state factorises with the hnal bath state being a 
thermal state, see [Appendix B Thus the bath B must remain in a thermal state under 
the optimal work extraction process. This clearly hints in the direction of the bound 
being reachable in the classical macroscopic limit as was suggested in [T]. 

Now one can consider special cases of this general relation. A subclass of map 
Eq. (36) are the single level to single level transitions, sketched in Fig. and in Fig. |^. 
The necessary condition for a transformation ps —)■ p's to be thermodynamically allowed 
here becomes 


(to) F(ps) - F(p's). 


(39) 


To maximise w one clearly wants to choose the hnal state thermal, p'g = ts, resulting 
in the highest permissible value = F{ps) — F{ts). In general one has Ff^^^{ps) < 
F{ps), and hence with approaching F{ps) in the i.i.d. limit 

(inhnitely many identical copies) [1]. This means that inequality Eq. (39) is here a 
lesser requirement than the tight bound Eq. ([^, giving just an upper threshold of what 
transformations may be allowed. 

For weight transitions from a single level to multiple levels, sketched in Fig. [^, the 
inequality gives 


(w) ^ F(a‘„) < F(ps) - F(p‘s). (40) 

So it is down to the choice of a suitable hnal weight state a'^ to enable the transformation 
to be thermodynamically possible. Fig. i) shows an example of a hnal energy 
distribution of the weight (orange box). It can be seen that if one was to shift up 
the hnal distribution in energy while keeping the entropy the same, F{a'y^) would grow 
thus making the inequality continuously tighter. The limit on how far the hnal weight 
distribution can be shifted in energy and the transition still being allowed is given by 
the system’s free energy diherence. 

For a multiple level to multiple level transition, sketched in Fig. [^, an example 
initial weight distribution is sketched in Fig. |^. If the energy distribution of the weight 
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is solely shifted to a higher energy without changing its shape or entropy, as indicated 
in blue in Fig. [^, then the inequality reduces to 

W £ < F(ps) - F(p's). (41) 


This result coincides with the tight bound derived for a more restrictive process in which 
global unitaries are chosen to be of a specihc form but for which only average global 
energy conservation is required na. 

Condition Eq. (38) produces sensible results for the special cases discussed here. 


While for the transformation to be thermodynamically possible fulhlling the condition 
is only necessary - not sufficient - the condition’s strength lies in its applicability for 
general globally energy conserving transformations on initially factorised states. Apart 
from providing an easy-to-check criterion to rule out many potential weight transitions, 
calling the transfer-quantity “work” can be supported, see e.g. [H]. Indeed we observe 
that condition Eq. (38) for the transfer quantity {w) is of the same form as the well- 


known bound on the average work (W) that can be extracted from a statistical ensemble 
when changing its probability distribution ps to another probability distribution p'g, e.g. 
in the context of Jarzynski’s equality [3], 


{W) <-AFs, 


(42) 


which is widely viewed as the second law of thermodynamics. 


8. Discussion 

We have re-derived the extractable work and the work of formation in the single 
shot setting proposed in [T], where a work storage system (a) transfers from a single 
energy level to another single energy level with probability 1 — e. Our alternative 
derivation of the optimal work value for which such a transformation is still possible 
was based on the comparison of subspace dimensions while no discussion of thermo- 
majorisation was required. (Thermo-majorisation pQ offers an additional and separate 
method of addressing the same optimisation problem.) The approach presented here 
facilitates a discussion of the hnal state of the system after work extraction has been 
performed. Indeed, we hnd that typically the reduced hnal state is a minimum free 
energy thermal state as has been suggested [T] . “Typically” is here to be understood in 
the very same sense in which it is often used in the context of thermalisation studies: 
it is overwhelmingly frequent w.r.t. Haar distributed unitaries [28l |29l 1301 El]- The 
implication of this Ending is two-fold: on one hand it should not come as a surprise if 
a “work-extracted” system is left in a thermal state, on the other hand the question 
of how far away from a thermal state a system may end up for a specihc, for instance 
physically motivated, unitary is a promising direction for future research. 

The presented approach opened the possibility to discuss a single shot work concept 
when transitions of the work storage system between (b) a single energy level and a range 
of final energy levels is permitted. Extensions to multiple levels are desirable to rehect 
experimental constraints - current and near-future experimental control of meso/macro- 
scopic systems does not allow to distinguish between energy levels that are spaced by 
less than an optical phonon. While the results for the multiple hnal level situation 
are mathematically sound, they raise interpretation issues when applied to the physical 
context. Besides showing an apparent violation of the second law that requires careful 
attention, an important conclusion emerges. The term “work extraction” intuitively 
suggests that an energy of a certain amount or higher is being made available but the 
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“minimum extractable work” found for case (b) was in fact higher than the single level 
case (a) and thus at odds with physical intuition. The single shot work analysis thus 
turns out to be not a good setting for discussing “work extraction”. It is however rather 
suited to characterise resonance processes where a specific amount of energy needs to be 
stored in the work storage system with high probability. Due to its nature, the energy 
stored may here be thought of as resonance work. 

Finally, the aim was to consider multiple energy levels of the work storage system 
for both, initial and final state. We have not tried to generalise the single shot work 
extraction approach to this situation because of the issues in physical interpretation that 
we could not fully resolve. Instead we provide a means to decide whether a transfer can 
be thermodynamically allowed or not by introducing a transfer-quantity (w) that we 
chose as the free energy of the weight. The derived necessary bound, {w) < —AFs, is 
applicable for transformations with general initial and final states of the weight. The 
bound is not tight in general, i.e. there are transformations that fulhl the bound while 
still not being possible with the thermodynamic resources considered here. However, 
the bound gives tight results for the single level case (a) when applied to iV —)■ oo 
independent identical (i.i.d) copies [H [15] and to a set of restricted global unitaries that 
results in multi-level to multi-level weight transfers na. Since the criterion is fairly 
simple to apply it facilitates the discussion of many practical limits on work extraction 
in the quantum regime. To rephrase, while Eq. (38) is not precisely tight it is sufficiently 
tight to be relevant and thus it may be considered a valuable concept. 
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Appendix A. Derivation of properties of the final states 


This part of the Appendix is dedicated to a more detailed derivation of the resnlts on 
the final states discnssed in Section HI 

We start by denoting the projector spanned by all eigenstates of the initial state t] 
that have total energy E and h{Es, g) = Ihy n®j(e, rf). The dimension of this projector 
dkffe^rj) = tr[nj^j(e,7])] has already been calcnlated in (14). By construction it equals 
the dimension d^^ of the respective projector n|(^. The maximum work extraction itself 
simply consists in a unitary V which maps the initial subspaces on the respective final 
ones, thereby respecting energy shells, which is the only condition on V. As mentioned in 
Sect. I^the crucial point here is that there is obviously not a single, unique unitary that 
corresponds to maximum work extraction but infinitely many. Since there are very many 
unitaries there may be potentially very many different reduced states of the system after 
maximum work extraction. In order to better understand the multitude of final system 
states it is instructive to consider unitaries which are randomly distributed according to 
the Haar measure. (Those unitaries may be constructed by drawing (column) vectors 
component-wise at random and doing a Gram-Schmidt orthonormalization For more 
details on averaging over Haar distributed unitaries see m- In the following we outline 
some results on distribution of quantities generated by such random unitaries without 
(formal) proof. 

Consider the unitary mapping V of some matrix A and a subsequent projection onto 
its diagonal part w.r.t. some orthonormal basis {|n)}, i.e., W){n\VAV^\n){n\. The 
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(n\ 


(A.l) 


average over all such unitaries may simply be inferred from symmetry considerations: 
{{^ \n){n\VAV^\n){n\)) = fl, fl := ^ \n){, 

n=l n=l 

where ({■)) denotes the average over all unitaries and d is the dimension of the 
Hilbertspace, d = tr[n]. Denoting the sum over m(< d) of the diagonal elements of 
this operator by s, we hnd for the latter: 


s : = 


J2HVAV^\n), 


n=l 


we find from (A.l) for the average of s: 

mtr[A] 


((^)) = 


d 


(A.2) 


(A.3) 


It is more subtle to compute the variance of s but it has been done, e.g., in the context 
of quantum typicality [29l [31]. From the latter works we find: 


(A^)^ := ((^^)) - ((^))^ = 


m(l - f)tr[A]^ 


(A.4) 


d{d+l) 

and for the relative deviation (setting d 1 ^ d for simplicity) 

Thus, if m is sufficiently large, those relative deviations become (negligibly) small. 
Which means even for a single unitary drawn at random according to the Haar measure 
s may be expected to be given by s ~ {{s)) for the overwhelming majority of cases. 

Equipped with these findings we may now identify the “typical” local reduced final 
state of the system S. To this end we consider (A.l) and make the identifications: 

/I :n®i(€,r,),,ng,(£,,,), n:nL (A.e) 

furthermore we note that 

tr[A] = tr[flA(e, r]) r] Ilgfe, t])] = (1 - e) PA. (A.7) 

Let be the diagonal part (w.r.t. the product energy eigenstates of the uncoupled 
subsystems) of the full system final state that has actually changed under work 
extraction, i.e. everything except for the e-part that remained unchanged, and ^ the 


state’s projection into the subspace of energy E. Then one finds from (|A.1|), inserting 

(1 


(|A.6|) and (|A.7|): 


( (^diagji?) ) 


ei P® TT^ 

-'ini -'-'-fin 
“fin 


(A.8) 


Since PF is invariant under V it may be calculated most conveniently from (15). From 
this calculation it is found to be actually independent of P, i.e., the total probabilities 

^ ^ is 

w) ® nvv('“’)- 


on all energy shells are the same. Using (|10|) and (|11|) yields 

( iVdiagjE/ 


TTe ^s{Es) ® nB(P - Es 
Es 


(A.9) 
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At this point it is actually more convenient to go back to the original, local energies and 
sum over £' 5 , Eb rather than Es, E. Doing so yields: 

((^kg)) « Us{Es)e-^^^ ® flsiEB - w)e-^^^ ® fl^(n;) (A.IO) 

Es,Eb 

Using the index shift Eb — w —?• Eb, this may be summed as 

((^diag)) 0^ (8 ® (A. 11) 

which is as given and discussed in Section]^ 

However, to conclude that the actual final state of the system S corresponding to 
a single random unitary is close to e~^^^ /Zs and thus a minimum free energy state it 
remains to be shown that the relative deviations of are small. To this end it is 
instructive to write out the diagonal matrix elements of the reduced final system state 
0-5 = tYBwW] explicitly 


{Es,g\o's\Es,g) = '^P^n{E,Es,g), (A. 12 ) 

E 

Mb{Eb) 

Pfin{E,Es,g) := (Es,g,E,f,w\r]'^i^g\Es,g,E,f,w), 

/=i 

where \Es, g, E, f,w) are the respective (product) energy eigenstates of the decoupled 
system. So a single Pa^^E, Es, g), i.e., a contribution to the respective diagonal element 
of the final reduced state of S from some energy shell E, is itself a sum over very many 


diagonal elements of the full final state, as may be seen from the lower line of (A. 12 ) 
Namely, it runs over all bath energy eigenstates at Eb = E — Es — w, i.e 
Mb{Eb) summands. Thus from (A.5), setting m = Mb{E 


contains 
it may be inferred 


Es — w 

that the relative deviations of the contributions to the diagonal matrix element are on 
the order of 

^PRn{E, Es, g) 1 


{{PUE,Es, 9)}) yMB(E-Es-w) 


(A. 13) 


which is very small if the bath is large. This is (22) as given and discussed in Section 
^ Furthermore, if the unitaries are drawn at random, the Ue are independent of each 


other. Thus one may conclude from (A. 13) and the central limit theorem, that according 


to the first line of (A. 12 ), the relative deviations of {Es, g\<Js\Es, g) will be even smaller 
if many energy shells are involved. 

An analogous consideration (which we do not present here in full detail for brevity) 
appli es to the off-diagonal elements of the final reduced state of S. The analogue to 
(A.3) required to find the average of those off-diagonal elements reads [ 21 ] 

{{{n\VAV\n'))) = 0 for n 7 ^ n'. (A.14) 


Thus, using the same identifications as before (A. 6 ), on average off-diagonal elements of 
the full system final state vanish. Since off-diagonal elements of reduced states are sums 
of off-diagonal elements of the full state exclusively, the averages of the final reduced 
state of the system as also vanish. Thus on average, the final state is indeed diagonal 
in its energy eigenbasis. Furthermore only those (off-diagonal) elements of the full state 
for which the full system states In'), |n) both correspond to the same bath system state 
\Eb, f) contribute to the reduced state of S at all. Thus each contribution to the reduced 
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state of S may be associated with a certain energy shell Eb oi the bath. An analogue 

dements 

for n 7 ^ n', (A. 15) 

(A.16) 


to (A.4) quantihes the deviations of off-diagonal elements, 

(Ao7 :=((|{n|V'/lV'tK>n) 


and reads, according to e.g. 

tr[A2] 


(Ao)^ : = 


d{d 1) 


From (A.6) and (A.l) it may be inferred that, tr[A^] < tr[fl;^j(e,r^)] and d = trpgj. In 
Sect. is has been established that both traces essentially scale with the dimension of 
the bath, i.e. 

trlng|(€, I,)] oc Mb{E), trlnfj k Mb{E). (A.17) 

Thus we find 

C 


Ao < 


(A. 18) 


\/Mb{E)' 

where C is a number which is independent of the overall dimension of the bath. As 
Mb{E) > Mb{E — Eg — w), this implies that the deviations of the off-diagonal elements 


of the reduced state of S again scale at most as 


= , just like in the case of 


Mb{E—Es—w) ‘ 

the diagonal elements, cf. (A.13). This result is also given and discussed in Section]^ 


Appendix B. Derivation of upper bound to the state transfer quantity (w) 


This part of the appendix is dedicated to the derivation of (38) in Sect. In the 
remainder any state of the total system, possibly including correlations, will be denoted 
by T]. Reduced local states will be denoted by ps = tiiBwiv] fhe non-equilibrium 
system, aw = tr^^)?]] for the weight storage system, and tb = sw[v] for the bath. 
Just as before V transforms the total initial state rj into the total final state rj' which 
means p' = VpV"^. And just as before we assume conservation of the sum of local 
energies, i.e., [Hg Hb + Hw, R] = 0. Free energies will be denoted as 

S{X) 


EiX) := U{X) - 


/S 


with S{X) := —tr[X InX], 


(B.l) 


where X := Ps,(^w,tb,P- Thus, free energies may refer either to parts or to the total 
system. The expectation value of the respective energy is denoted by U{X) for the 
corresponding Hamiltonian, Hg, Hw, Hb- Generally primed operators will refer to final 
states and unprimed operators to initial states. 

Again, as also mentioned in Sect. we assume factorizing initial conditions, i.e., 
p = Ps ® ® CFw- We also make the assumption of the initial state of the bath is 

thermal, tb oc The only difference in the set up of the approach in this Section 

compared to Sections [3]j^ is that we do not require the state density of the spectrum of 
the bath to be exponentially growing. By construction, the total entropy and the sum 
of local energies are both invariant under V. Thus, the initial total free energy and the 
final total free energy are the same, E{p) = E{p'), and due to the factorizing initial 
conditions we have 


E{p) = E{pg) + E{tb) + E{aw)- 


From the Araki-Lieb theorem [IB] it follows that 


S{p') < Sip's) + Sir'B) + Sia'^) 


(B.2) 

(B.3) 
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and hence, 


F{^')>F{p's) + F{T'^) + F{a'^). 


Combining (B.2) and (B.4) yields 


(B.4) 


(B.5) 


(B.6) 


F(ps) + F{tb) + F{aw) > F{p's) + F(t'b) + F(a'^.) 

which may simply be rearranged as 

F{a'^) - F{aw) < -{F{p's) - F{ps)) - {F{t'^) - F{tb)). 

A thermal state has the lowest free energy given a specihc inverse temperature fd. 
Since we assumed an initial thermal state for the bath it follows that F{tb) is the lowest 
possible bath free energy. Consequently, whatever F{t'q) is, one has 

- [F(Fb) - F(tb)) < 0, (B,7) 


By dropping the corresponding term in (B.6) we only make the r.h.s. possibly larger 


Since this is not in conflict with the inequality in (B.7) we may simply drop the term, 
obtaining 


FW^,) - F(a„) < -(F(p's) - F(ps)). 


(B,8) 


The l.h.s. is F{a'y^) —F{aw) = {w) as dehned in (37) completing the derivation of (38). 








